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Abstract 

We look for differential equations of the form 



where the coefficients {ci(x)}°^ are continuous functions on the real line and where 
{ci(x)}°l 1 are independent of n, for the generalized Jacobi polynomials {-P"^' M,Ar (aO}^ 
and for generalized Laguerre polynomials {L"' M > N (x)}^^ which are orthogonal with 
respect to an inner product of Sobolev type. 

We use a method involving computeralgebra packages like Maple and Mathematica 
and we will give some preliminary results. 
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1 Introduction 



In this paper we consider the polynomials \ P®'P' M ' N (x)\ (see 0) which are orthogonal 

I. J n=0 

on the interval [—1,1] with respect to the weight function 
!> + /? + 2) n , p 



2 «+/3+i r ( a + !) r ^ + x) 



1 + xf + M5(x + 1) + AT«5(a; - 1), 



where a > -I, (3 > -1, Af > and iV > 0. 

We also consider the polynomials \ L®' M ' N (x) \ (see ||) which are orthogonal with 

I J n=0 

respect to the Sobolev inner product 

oo 

<f,9> = r , \ u I x a e- x f(x)g(x)dx + Mf(0)g(0)+Nf'(0)g'(0), 
1 (a + lj J 

where a > -1, M > and iV > 0. 

We are looking for differential equations of the form 



oo 
i=0 



Ci(x)yW(x) = 0, 



satisfied by these sets of orthogonal polynomials, where the coefficients {ci(x)}°l are contin- 
uous functions on the real line and {ci(x)} c *L 1 are independent of the degree n. 

In [^j J. Koekoek and R. Koekoek found a differential equation of the above type for the 

generalized Laguerre polynomials \L" ,M (x) \ (see [f7|) for all a > —1. These polynomials 

l. > n=0 

form a special case of the above mentioned Sobolev Laguerre polynomials, since L"' M '°(x) = 
L^' M (x). For more details the reader is referred to [|j and section 2 of this paper for more 
results concerning the coefficients of this differential equation. 

The differential equation found in j|] was computed by hand, without the help of comput- 
ers. This is nearly impossible for these other cases. We need computers to handle the very 
huge expressions we have to deal with. 

The method we used to find the results listed in this paper is explained in [|]]. We refer 
to this report for more details. 

2 The infinite order Laguerre differential equation 

In 0| the following theorem was proved without the use of computers : 

Theorem 1. For M > the polynomials \ L?!' M (x)\ satisfy a unique differential 

I J 71=0 

equation of the form 

oo 

M^2ai(x)y {i) (x) +xy"{x) + (a + 1 - x)y'(x) + ny(x) = 0, (1) 

i=0 

where {cii(x)}°^ are continuous functions on the real line and {aj(x)}^ 1 are independent of 
n. 
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Moreover, the functions {ai(x)}°^ are polynomials given by 

Oq(x) 



n + a + 1 
n — 1 



(2) 



= I (J * J J (j^j (« + 3)i-M < = 1, 2, 3, . . . . 



Later we discovered that the coefficients have the following interesting property. 

Theorem 2. The coefficients {ai{x)} c *i l of the differential equation given by and (^) 
satisfy 

Sin7TQ X „ / 1 



Ti 



.. vr (a + 2)(a + 3) ^ a + 4 

For nonnegative integer values of a we have : 



x , a > —1. 



£(*) a<(x) = (-l) Q+fc a fc (-x), fc = 1,2,3,.... 



i=k 



Note that this theorem implies for nonnegative integer values of a : 

oo oo 

y^aj(aj) = and ^zaj(x) = (— 
i=l i=l 

The proof of this theorem can be found in GJ . 

3 Some preliminary results for the Sobolev Laguerre polyno- 
mials 

In this section we look for a differential equation of the form 

oo oo 

M ]T Oj {x)y® (x) + NJ2 k{x)y {{) (s) + 
i=0 i=0 

oo 

+ MN^2ci(x)y {i) (x) +xy"(x) + (a + 1 - x)y'(x) + m/(x) = (3) 

i=0 

for the polynomials 

y(x) = L% M > N (x) = A L^\x)+A 1 ^-L^(x)+A 2 -^Li a \x), 
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where the coefficients Aq, A\ and A<i are defined by 
A 



i+M n+a ) + n( " +2 \: (a + 1) jv( n+a ) + 

Vn-ly (a + l)(a + 3) \n-2J 



+ 



MN n + a\ n + a + 1 

(a + l)(a + 2)L-lM n-2 



At = M 



n + a 



n 



(n-1) fn + a\ 2MN (n + a\/n + a + V 
(a + 1) In - 1 J + (a + 1) 2 I n Jl n-2 



(4) 



A 2 



iV /n + a 



+ 



MAT /n + a\ /n + a + 1 



(a + l)U-li (a + 1 



n — 1 



For details concerning these generalized Laguerre polynomials and their definition the reader 
is referred to || and |J. 

Of course, since L" ,M,0 (x) = L® ,M (x), the coefficients {aj(x)}°^ are given by (|2j). 

Although the general form is still an open problem so far, we know that the differential 
equation given by (^) is not unique as in the case of the differential equation ([!]). This is 
explained by the following theorem : 

Theorem 3. The polynomials \ L" ,M ' N (x)X satisfy the following infinite order differ- 

L J 71=1 

ential equation : 

oo oo 

£ b* (a, x)y® (x) + M £ <(a, x)y« (x) = 0, 



i=0 



i=0 



where 



and 



b*(a,x) 



*-3 



x j , i = 0,1,2,... 



c*(a,x) 



-1)*' 



x 4 , i = 0,1,2,... . 



The proof is very easy and is based on the observation that 

-n) k 



^6*(a,x) J D^4«)(x) 



i=0 



nT(k) 



n > 1, fc = 0,1,2 



and 



fc*(a,x)^4«)(x)= ( n + a )±^-, /,-(). 1.2. 



i=0 



n /(a + 1 



By using these formulas and the definition (||) of the coefficients Aq, A\ and A^ the result 
follows immediately. 

Again the reader is referred to the report Q for more details. 

In the special cases a = 0, a = 1 and a = 2 differential equations of the form dH) are 
found explicitly. In these three special cases of integer values of the parameter a we find a 
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linear differential equation of formal order 4a + 10. By formal order we mean that for special 
cases (M = or N = 0) the true order might be lower. We will not give all results here, but 
we again refer to J|] for all details. 

As an example we give the result in the case that a = 0. In that case we have found the 
following linear differential equation of formal order 10 : 

-^MiVx 5 y (10) (x) + ^MN{5x A - x 5 )y (9) (x) + 



+ 
+ 
+ 
+ 
+ 

+ 



— MN(72x 3 - 45x 4 + 4x 5 ) - —Nx 4 
24 v 7 12 



y (8) (x) + 



y {1) {x) + 



^MiV(36x 2 - 72x 3 + 20x 4 - x 5 ) - ^N(4x 3 - x A ) 

^MN(-252x 2 + 224x 3 - 35x 4 + x 5 ) + ^iV(-10x 2 + 9x 3 - x 4 ) 
1 1 

— MN(1680x 2 - 840x 3 + 75x 4 - x 5 ) + -iV(-12x + 81x 2 - 33x 3 + 2x 4n 
60 6 



y^(x) + 



y (5) (x) + 



1 

24 



M7V(-420x 2 + 120x 3 - 5x 4 ) + 



+ — iV(24 + 36x - 150x 2 + 34x 3 - x 4 ) - -Mx 2 
12 v 1 2 

-M7V(15x 2 - 2x 3 ) + -N(-6 - 2x + 9x 2 - x 3 ) + M(-2x + x 2 ' 
3 2 

111 

--MNx 2 + -N(2 - x 2 ) + -M(6x - x 2 ) + x 



y (4) (x) + 



y (3) (x) + 



y"{x) + [1 — (M + l)x] y'{x) + 



H n 

120 



MN(n 2 - l)(n + 2)(2n + 1) + 10iVn(n 2 - 1) + 60M(ra + 1) + 120 y{x) = 0, 



for the polynomials 



y(x) := L% M > N (x) = A L n (x) + A x L' n {x) + A 2 U' n {x) 



with 



and 



L n (x) :=4°)(x) = £ 



k=0 



(-l) k (n> 
~k\ i k , 



x k , n = 0, 1,2, . . . 



1 1 

A = 1 + Mn + -Nn(n - l)(2n - 1) + —MNn 2 (n 2 - 1) 
6 12 

Ai = M + Nn(n - 1) + -MNn(n 2 - 1) 
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A 2 = Nn + -MNn(n + 1). 
The results obtained in these three special cases give rise to the following conjecture. 



Conjecture. If a is a nonnegative integer, the polynomials { L^ ,M ' N (x)\ satisfy a 

I J n=0 

differential equation of formal order 4a + 10 which is of the form 



2q+4 2a+8 

M £ a l (x)y^(x)+N ]T 0i(x)y®{x) + 
i=0 i=0 
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4q+10 

+ MN Ji( x )y^( x ) + xy"(x) + (a + 1 - x)y'{x) + ny{x) = 0, 

i=0 

where the coefficients {^(x)} 2 ^ 4 , {/^(x)} 2 "^ 8 and {li{x)} 4 =i 10 are polynomials independent 
of n which satisfy 

2o+4 2o+8 4o+10 

i=l j=l i=l 



4 The generalized Jacobi polynomials 

In this section we will deal with the problem of finding a differential equation for the gener- 
alized Jacobi polynomials \P"'^' M ' N (x)\ 

Since the well-known second order differential equation for the classical Jacobi polynomials 

| p( a 'Ph x \\ i s given by 

I J n=0 

_2\ /// 



(1 - x 2 )y"{x) + [/3 -a- {a + /3 + 2)x] y'{x) + n{n + a + (3 + l)y{x) = 0, 
it is clear that we look for a differential equation of the form 

oo oo oo 

M^ai (x)y^ (x) + NY k (x)y (t) (x) + MN £ a (x)y^ (x) + 

i=0 i=0 i=0 

+ (1 - x 2 )y"{x) + [(3 -a - {a + (3 + 2)x] y'(x) +n(n + a + (3 + l)y(x) = 0, 

where 



and 



y {x) = P^ M > N ( X ) = A P^\x) + [Ai(l - x) - A 2 {1 + x)] ^^(l) 



(n + (3\(n + a + f3+l\ (n + a\(n + a + (3 + l\ 
In — 1/1 n I In— 1/1 n / 

A = l + M— J -\ L + N \ L\ L + 



+ MN 



(a + P + 2) 2 (n + a + [3+l\ 2 



{a + 1)(J3 + lH n-1 



n + /3\/n + a + /3 



^ _ M \ n )\ n J MN (n + a + fi\(n + a + f3 + l 
1 ~ {a + (3+ I) rn + a\j + JoT+T) [ n-1 ) \ n , 

(n + a\(n + a + (3\ 
N I, " j\ n J MN (n + a + f3\fn + a + /? + 1 N 
2 " (a + 0+1) Tn~+~p\ + W+l){ n-1 ){ n , 



Here we used the same definition as in Q, but in a slightly different notation. 

In this case the differential equation will be unique in its general form if it exists. We 
introduce the notation 



ao(x) 

< 


= ao(n, a, (3, x) 


ai(x) 


= ai(a,f3,x), i 


[ b (x] 


= b (n,a,/3, x) 


\ h{x) 


= bi(a,/3,x), i 


{ c (x) 


= c (n,a,f3,x) 


{ Ci(x) 


= Ci(a,f3,x), i -- 



and 



Since the polynomials \ P£>P> M ' N (x)\ satisfy the symmetry relation 

I J n=0 

pa,P,M,N(_ x) = (_ l)np P,a,N,M( x) 

we have 

ao(n, a, /3, — x) = bo(n, ft, a, x) 
a,i(a, f3, -x) = (-lfbiiP, x )> i = 1, 2, 3, . . . 
co(n, a, (3, —x) = co(n, 0, a, x) 



and 



Ci(a, (3, -x) = (-l)Vj(/3, a,x), i = 1, 2, 3, 

The general form is still an open problem, but the symmetric case (3 = a and N = M 
turns out to be much less difficult. If f3 = a and N = M we can choose Cj(x) = for all 
i = 0,1,2,.... 

In fact we have shown that there exists a differential equation of the form 



M^ ai {x)y^{x) + (1 -x 2 )y"{x) - 2(a + l)xy'(x) + n(n + 2a + l)y(x) = 0, 



i=Q 



where 



and 



//!./•; - P n tt ' Q ' M ' M (x) = C Q P^ a \x) - C lX ^-P^' a \x) 



d 



(5) 
(6) 



Cn = 1 + M- 



2n fn + 2a + l\ r9 ( n + 2a + 1 



a + 1) 



+ AM Z 



n 



n — 1 



(7) 



Ci 



2M n + 2a\ 2M 2 ( n + 2a \ in + 2a + 1 



+ 



n 



(a + 1) V n - 1 



(2a + 1) 

This differential equation turns out not to be unique. We write 

ao(x) := ao(n, a, x), n = 0, 1, 2, . . . 

a>i{x) := aj(a, sc), i = 1, 2, 3, . . . . 



n 
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If we substitute (jfjj) and (|?|) in the differential equation then we finally find three 
equations for the coefficients {aj(x)}°^ which are equivalent to the following two : 



J^a i (x)D i P^ a ' a \x) 



i=0 



4 (n + 2a\ d 2 ^ 
(2a + 1)1 n dx 2 n K } 



and 

£ m,(x)Z^)(x) + xf a % {x)D^Pt a \x) = J n + + ^ ^P^H*)- 

i=0 i=0 \ / 

We introduce the notation 

ao(n, a, x) = ao(l, a, x)&o(?"i<> a, + co(n, a, x), n = 0, 1, 2, . . . 

aj(a, x) = ao(l, a, x)6j(a, x) + Cj(a, x), i = 1,2,3, ... . 

Now ao(l, a,x) is arbitrary, but {&i(x)}°^ and {cj(x)}^ are uniquely determined. The 
explicit form of these coefficients is given in the two theorems below. 



Theorem 4. The polynomials \ P^ ,a ' M,M (x)\ satisfy the linear infinite order differ- 

L J n=0 

ential equation given by : 

oo 

5>(x)y«(x)=0, 



i=0 



where 



6 (x) := b (n, a,x) = -[l- (-l) n ] , n = 0, 1, 2, 



51-1 



6i(x) := 6j(a,x) 



-(-x) i , i = 1,2,3,.... 



Now we come to our main result 



by: 



Theorem 5. The polynomials { P^ ,a ' M,M (x)\ satisfy the differential equation given 

y- > n=0 



M ]T Ci(x)y W (x) + (1 - x 2 )y"(x) - 2(a + l)xy'(x) + n(n + 2a + l)y(x) = 



where the coefficients {ci(x)}°l Q are defined by 

c (x) := c (n, a, x) = 4(2a + 3) ( n 2 J , n = 0, 1, 2, 



and 



Ci (x) = (2a + 3)(1 - x 2 )c*(x), i = 1, 2,3, ... , 
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where 



c\{x) := c\(a, x) = 



4i x ) ■= c*i{a,x) 



* = 2,3,4, 



I" 1 n 



We remark that this theorem implies that if a is a nonnegative integer, the polynomials 



a,a,M,M i 



n=0 



satisfy a linear differential equation of formal order 2a + 4. 



Finally, we note that since 



and 



we have 



2-^1 



2-^1 



—i + l,a + 



L| ^ (-°-2 + 0«-i >i = 1>2>3> _ 



-i + 1, a + 



(5)*-! 

(-a - 1 + 
(§)i-i 



, * = 1,2,3,... 



E4(«--i) = E71 



(-a - (-a - 2 + f 



i-l 



i=l 



i(i)i(l)i4< (i-l)! (i)<-i 



<I(|)i(2)i(i)i 



^ (2*)»(3)2i 



and 



E c 2i+i( a ' _1 ) 



E 



i=0 



(a + 1) (-a)i-i (-a - 1 + 
t (l)i(|)i4< (i-l)! (§)<_! 

_l( a + 1) J ( -a )i( -a + *)* , 1 u 



<l(2)i(§)i(|)i ^ 



(« + i)E 



(-u) 2 i n2 i 



= 2^ (~ a ~ iWl 2 2i+l 



-(2i + l)!(4) M ^(2i + l)!(3) 2m 



In the same way we find 



OO OO / -i \ 



i=l 



^ (2»)K3)a 



and 



This implies that 



OO OO / -i \ 

E 4 + i(«, i) = -2E ^ ( ~ ~ 2m 22m - 



i=0 



-(2z + l)!(3) 2m 



^c*(a,-l) = ^4(a,-l) + E c 2i+i («>-!) 

i=l i=l i=0 

_ ^ (-«-l)i i_o p f 
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and 



oo oo oo 

8=1 8=1 8=0 



i=l 
oo 



i=0 



i!(3)i 



-a — 1 
3 



2 . 



Remark. The proof of theorem 5 will be given in a forthcoming paper ||. Further we 
remark that this result proves one of the conjectures by L.L. Littlejohn and W.N. Everitt 
given in HI. 
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